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Abstract 

Deterministic equilibrium flows in transport networks can be investi- 
gated by means of Markov's processes defined on the dual graph represen- 
tations of the network. Sustained movement patterns are generated by a 
subset of automorphisms of the graph spanning the spatial network of a 
city naturally interpreted as random walks. Random walks assign abso- 
lute scores to all nodes of a graph and embed space syntax into Euclidean 
space. 
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1 Motivation 

In most of researches devoted to the improving of transport routes, the opti- 
mization of power grids, and the pedestrian movement surveys the relationships 
between certain components of the urban texture are often measured along 
streets and routes considered as edges of a planar graph, while the traffic end 
points and street junctions are treated as nodes. Such a primary graph repre- 
sentation of urban networks is grounded on relations between junctions through 
the segments of streets. A usual city map based on Euclidean geometry can be 
considered as an example of primary city graphs. 

On the contrary, in another graph representation of urban networks known 
as dual, relations between streets are encoded through their junctions. Streets 
within that are treated as nodes while their intersections are considered as edges. 



Dual city graphs are extensively investigated within the concept of space syntax, 
a theory developed in the late 1970s, that seeks to reveal the mutual effects of 
complex spatial urban networks on society and vice versa, [2 [5] . 

Recently, while acting as a referee in Physica A, we have been asked to review 
a manuscript devoted to city space syntax in which the dual graph representa- 
tion of urban patterns is used. Among the technical remarks expressed by the 
previous referees, we have discovered the following criticism: 

"[The] dual option in representing cities would have deserved a 
deeper discussion, as it must be said that no one in the world, space 
syntax apart, appears to use it ... space syntax supporters and 
scientists should dedicate more efforts to the legitimation of their 
theory's fundamental notions in the universally accepted language 
of the wider scientific community. " 

Being not obliged to reply personally, we nevertheless must answer such a 
corporative challenge by explaining the deep arguments in favor of studying 
namely the dual graph representations. 

Legitimacy of the space syntax approach is based on the surprisingly high 
degree of correlation between aggregate human movement rates and spatial con- 
figuration analyzed in space syntax theory by means of dual city graphs. We 
recall here that the key result of space syntax research is that the pattern of spa- 
tial integration in the urban grid is a key determinant of pedestrian movement 
in cities across the world, 3j. 

We emphasize that the dual graph analysis neither incorporates many of the 
factors considered important in previous models of human behavior in urban 
environments (such as the motivations and the origin-destination information) 
nor direct account was taken of the metric properties of space (lengthes of 
streets), [3]. Nevertheless, the robustness of agreement between space syntax 
predictions and rush hour movement rates is now supported by a number of 
similar studies of pedestrian movement in different parts of the world and in 
an everyday commercial work of the Space Syntax Ltd., [5]. Similar results also 
exist for vehicular movement 6] showing that the spatial configuration of the 
urban grid is by itself a consistent factor driving transport flows. 

The high quality agreement with the empirical data obtained from pedestrian 
movement surveys looks amazing even for those who invented the space syntax 
approach. Quoting ^jj there is nothing else to evoke but that 

"the way people understand their environment and decide on move- 
ment behaviors is somehow implicitly embedded in space syntax 
analysis". 

In the next section, we address three questions of essential importance for the 
analysis of transport networks: 

1. Why does the space syntax approach based on the pure configurational 
analysis of dual graphs appear to be so successful in prediction of pedes- 
trian and vehicular flows, essentially in rush hours? 
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2. Why does the measurement of the simultaneous importance of spaces in 
a city network by a quantity which is nothing else, but an element of 
a transition probability matrix determining a Markov chain fit best the 
empirical data on crowding patterns? 

3. What is the geometrical structure beyond space syntax? 

We show that equilibrium transport flows through edges of a connected 
undirected graph G at equilibrium are related to the stationary distribution of 
random walks defined on its dual counterpart G* . Random walks constitute a 
fundamental process defined on a graph generated by the set of its automor- 
phisms preserving the notion of connectivity which assign the absolute scores to 
all nodes of the graphs. Random walks embed the space syntax of any transport 
network into the Euclidean space E^^^ in which distances and angles get clear 
statistical interpretations. 

2 Mathematics of transport networks 

Any graph representation naturally arises as the outcome of a categorization, 
when we abstract a real world system by eliminating all but one of its features 
and by grouping together things (or places) sharing a common attribute. For 
instance, the common attribute of all open spaces in city space syntax is that 
we can move through them. All elements called nodes that fall into one and 
the same group V are considered as essentially identical; permutations of them 
within the group are of no consequence. The symmetric group Sat consisting of 
all permutations of N elements [N being the cardinality of the set V) constitute 
therefore the symmetry group of V . If we denote by C VxV the set of ordered 
pairs of nodes called edges, then a graph is a map G{V, E) : E ^ K CM.^ (we 
suppose that the graph has no multiple edges). If two nodes are adjacent, 
€ E we write i ~ j. 

2.1 Why dual graphs? 

First, we establish a connection between transport flows on the graph G and 
random walks on its dual counterpart G*. 

Given a connected undirected graph G{V, E), in which V is the set of nodes 
and E is the set of edges, we introduce the traffic function f : E ^ {0,oo[ 
through every edge e G E. It then follows from the Perron-Frobenius theorem 
[7] that the linear equation 



where the sum is taken over all edges e' € E which have a common node with e, 
has a unique positive solution /(e) > 0, for every edge e £ E, for a fixed positive 
constant ft. > and a chosen set of positive metric length distances t{e) > 0. 
This solution is naturally identified with the traffic equilibrium state of the 




(1) 
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transport network defined on G, in which the permeability of edges depends 
upon their lengths. The parameter h is called the volume entropy of the graph 
G, while the volume of G is defined as the sum 



Vol(G) ^\Y. ^(e) 



2 
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The volume entropy h is defined to be the exponential growth of the balls in a 
universal covering tree of G with the lifted metric, [5]-|ll|. 

The degree of a node i €V \s the number of its neighbors in G, degQ(i) = 
ki. It has been shown in that among all undirected connected graphs of 
normalized volume, Vo^G) = 1, which are not cycles and for which fci 7^ 1 for all 
nodes, the minimal value of the volume entropy, min(/i) = i X^iev ^'^S (^i ~ 1) 
is attained for the length distances 

. (2) 

2 mm(n.) 

where ki and kj are the degrees of the nodes linked by e e It is then obvious 
that substituting ([2]) and min(/i) into ^ the operator exp {—h£{e' j) is given by 
a symmetric Markov transition operator, 

/(e) = T , ^^''^ (3) 

where i and j are the nodes linked by e' G E, and the sum in ^ is taken over 
all edges e' £ E which share a node with e. The symmetric operator ^ rather 
describes time reversible random walks over edges than over nodes. In other 
words, we are invited to consider random walks described by the symmetric 
operator defined on the dual graph G*. 

The Markov process ([3]) represents the conservation of the traffic volume 
through the transport network, while other solutions of ([T]) are related to the 
possible termination of travels along edges. If we denote the number of neighbor 
edges the edge e £ E has in the dual graph G* as qe — degg* (e), then the simple 
substitution shows that w{e) = ^Jql defines an eigenvector of the symmetric 
Markov transition operator defined over the edges E with eigenvalue 1. This 
eigenvector is positive and being properly normalized determines the relative 
traffic volume through e G i? at equilibrium. 

Eq.Q shows the essential role Markov's chains defined on edges play in 
equilibrium traffic modelling and emphasizes that the degrees of nodes are a 
key determinant of the transport networks properties. 

The notion of traffic equilibrium had been introduced by J.G. Wardrop in [12j 
and then generalized in [13j to a fundamental concept of network equilibrium. 
Wardrop's traffic equilibrium is strongly tied to the human apprehension of space 
since it is required that all travellers have enough knowledge of the transport 
network they use. The human perception of places is not an entirely Euclidean 
one, but are rather related to the perceiving of the vista spaces (viewable spaces 
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of streets and squares) as single units and to the understanding of the topological 
relationships between these vista spaces, |14) . 

The use of Eq.Q also helps to clarify the inconsistency of the traditional 
axial technique widely implemented in space syntax theory. Lines of sight are 
oversensitive to small deformations of the grid, which leads to noticeably differ- 
ent axial graphs for systems that should have similar configuration properties. 
Long straight paths, represented by single lines, appear to be overvalued com- 
pared to curved or sinuous paths as they are broken into a number of axial lines 
that creates an artificial differentiation between straight and curved or sinuous 
paths that have the same importance in the system [TS]. Eq.Q shows that 
the nodes of a dual graph representing the open spaces in the spatial network 
of an urban environment should have an individual meaning being an entity 
characterized by the certain traffic volume capacity. 

Decomposition of city space into a complete set of intersecting open spaces 
characterized by the traffic volume capacities produces a spatial network which 
we call the dual graph representation of a city. 

2.2 Why random walks? 

While analyzing a graph, whether it is primary or dual, we assign the ab- 
solute scores to all nodes based on their properties with respect to a transport 
process defined on that. Indeed, the nodes of G{V, E) can be weighted with re- 
spect to some measure m = X^iey ^ij^ specified by a set of positive numbers 
TTLi > 0. The space £'^{m) of square-summable functions with respect to the 
measure to is a Hilbert space Ti-iV). 

Among all measures which can be defined on V, the set of normalized mea- 
sures (or densities), 



are of essential interest since they express the conservation of a quantity, and 
therefore may be relevant to a physical process. 

The fundamental physical process defined on the graph is generated by the 
subset of its automorphisms preserving the notion of connectivity of nodes. An 
automorphism is a mapping of the object to itself preserving all of its structure. 
The set of all automorphisms of a graph forms a group, called the automorphism 
group. For each graph G(V, E), there exists a unique, up to permutations of rows 
and columns, adjacency matrix A, the N x N matrix defined by Aij = 1 if « ^ j, 
and Aij = otherwise. As usual A is identified with a linear endomorphism of 
Co{G), the vector space of all functions from V into R. The degree of a node 
i £ V is therefore equal to 



Let us consider the set of all linear transformations defined on the adjacency 




(4) 




(5) 
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matrix, 

= ^^JsiAsi, T.jsi e M, (6) 

s,l=l 

generated by the subset of automorphisms of the graph G. 

The graph automorphisms are specified by the symmetric group Sat including 
all admissible permutations p € Sjv taking i e V lop{i) G V. The representation 
of Sat consists of all iV x matrices Hp, such that (Hp)^ ^j- .^ = 1, and (Hp) . ^ = 
if j ^p{i). 

The function Z (A)^^ should satisfy 

n;z(A) Hp = z(njAnp), (7) 

for any p € Sn, and therefore entries of the tensor must have the following 
symmetry property, 

•^P{i)p(j)p(s)p{l) ~ ^ijsli (8) 

for any p S Sat. Since the action of the symmetric group Sat preserves the 
conjugate classes of index partition structures, any appropriate tensor !F sat- 
isfying ([8]) can be expressed as a linear combination of the following tensors: 

{1, Sij,Sis,Sii,Sjs, Sji,Ssi, SijSjs, SjsSsi,SsiSii, SiiSij, 

Sij6si,6isSji,SiiSjs,SijSiiSis} ■ Given a simple, undirected graph G such that 
An = for any i £ V then by substituting the above tensors into ([6]) and 
taking account on symmetries we conclude that any arbitrary linear permuta- 
tion invariant function must be of the form 



Z {A)^- = ai+6ij {a2 + askj) + a4Aij, (9) 

with kj = deggQ) and ai.2,3.4 arbitrary constants. 

If we require that the linear function Z preserves the notion of connectivity, 

fc, ^ 5] Z(A)^^., (10) 

it is clear that we should take ai = 02 = (indeed, the contributions aiN and 
02 are incompatible with pO)) ) and then obtain the relation for the remaining 
constants, 1 — 03 = 04. Introducing the new parameter /3 = 04 > 0, we write 
(0) as follows, 

Z{A)^^ = {1-P)d,jk,+PA,,. (11) 

If we express (|10p in the form of the probability conservation relation, then the 
function Z (A) acquires a probabilistic interpretation. Substituting (|TT|) back 
into (fTO|l . we obtain 

= E,-ev(l-/^)'5..+/3 4^ (12) 
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The operator T^j represents the generahzed random walk transition operator 
if < /3 < k'^l^^ where fcmax is the maximal node degree in the graph G. In the 
random walks defined by , a random walker stays in the initial vertex with 
probability 1 — /3, while it moves to another node randomly chosen among its 
nearest neighbors with probability P/h. If we take f3 — 1, then the operator 
Tj;P describes the usual random walks discussed extensively in the classical 
surveys [TB]-[T5]. 

Being defined on a connected aperiodic graph, the matrix T^p is a real 
positive stochastic matrix, and therefore, in accordance to the Perron-Frobenius 
theorem [7], its maximal eigenvalue is 1, and it is simple. A left eigenvector 

TtT'-'^^ = TT (13) 

associated with the eigenvalue 1 has positive entries satisfying It is inter- 
preted as a unique equilibrium state tt (stationary distribution of the random 
walk). For any density a E Ti.{V), 

TT = lim crT*. (14) 

t — >oo 



2.3 Space syntax as Euclidean space 

Markov's operators on Hilbert space appear therefore as the natural lan- 
guage of complex network theory and space syntax theory in particular. Now 
we demonstrate that random walks embed connected undirected graphs into 
Euclidean space, in which distances and angles acquire the clear statistical in- 
terpretations. 

The Markov operator T is self- adjoint with respect to the normalized mea- 
sure ([4|) associated to the stationary distribution of random walks tt, 

f = i (^1/2 T ^-1/2 + ^-1/2 tT ^1/2^ ^ (15) 

where is the transposed operator. 

In the theory of random walks defined on graphs [16l [18] and in spectral 
graph theory [TO], basic properties of graphs are studied in connection with 
the eigenvalues and eigenvectors of self-adjoint operators defined on them. The 
orthonormal ordered set of real eigenvectors if^i, i — 1 . . . N, of the symmetric 
operator T defines a basis in 7i(V^). 

In particular, the symmetric transition operator T of the random walk de- 
fined on connected undirected graphs is 



T,, ^ y^' . (16) 




Its first eigenvector ipi belonging to the largest eigenvalue fii = 1, 

V'lf = Vl, i'h = ^^, (17) 
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describes the local property of nodes (connectivity), since the stationary distri- 
bution of random walks is 

TT,. = A (18) 



2M 



where 2M — J2iev '^^^ remaining eigenvectors, { ips }s=2^ belonging to the 
eigenvalues 1 > fi2 > ■ ■ ■ fJ'N ^ —1 describe the global connectedness of the 
graph. For example, the eigenvector corresponding to the second eigenvalue /i2 
is used in spectral bisection of graphs. It is called the Fiedler vector if related 
to the Laplacian matrix of a graph JT9]. 

Markov's symmetric transition operator T defines a projection of any density 
a £ Ti-lV) on the eigenvector -01 of the stationary distribution tt, 



(19) 



in which cr^ is the vector belonging to the orthogonal complement of ipi. In 
space syntax, we are interested in a comparison between the densities with 
respect to random walks defined on the graph G. Since all components ipi^i > 0, 
it is convenient to rescale the density cr by dividing its components by the 
components of "01, 

a.,= ^. (20) 



Thus, it is clear that any two rescaled densities a, p € H differ with respect to 
random walks only by their dynamical components. 



for all t > Q. Therefore, we can define the distance 
densities established by random walks by 



l|2 

pWt 



E 

t > 



ai_ - Pa 



T between any two 
(21) 



or, using the spectral representation of T, 



k - pWt Et >o ^J=2ys (o-^ - p-^\iJj){i'A^^ - p^) 

_ ( g-j_ -p^ \ i>s )( -03 1 -pj_ ) 



(22) 



i-p.. 



where we have used Diracs bra-ket notations especially convenient for working 
with inner products and rank-one operators in Hilbert space. 

If we introduce a new inner product for densities cr, p £ 7i(V^) by 



N 



>>P)l 



then ([^^ is nothing else but 



EE 

t>0 s=2 



[a±_ \ 'tl!s){'ips\p±) 
1 - A^s 



lkllT + ll/'llT^2(cr,p)j,, 



(23) 



(24) 
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where 

^^.f^^f-i^immM. (25) 

is the square of the norm of cr G Ti-lV) with respect to random walks defined 
on the graph G. 

We finish the description of the {N — l)-dimensional Euchdean space struc- 
ture of G induced by random walks by mentioning that given two densities 
(T, (0 G Ti.{V), the angle between them can be introduced in the standard way, 

cosZ(p,a) = ^fl^l^. (26) 

Random walks embed connected undirected graphs into the Euclidean space 
]g7V-i ^j^jg embedding can be used in order to compare nodes and to construct 
the optimal coarse-graining representations. 

Namely, in accordance to l|25|) . the density 6i, which equals 1 a.t i & V and 
zero otherwise, acquires the norm || Si Wj, associated to random walks defined on 
G. In the theory of random walks [16], its square, 

11^.11^ = ^Et^' (27) 

gets a clear probabilistic interpretation expressing the access time to a target 
node quantifying the expected number of steps required for a random walker to 
reach the node i G V starting from an arbitrary node chosen randomly among 
all other nodes with respect to the stationary distribution tt. 

The Euclidean distance between any two nodes of the graph G calculated as 
the distance ((22|) between the densities Si and Sj induced by random walks, 

Kr,j = \\d,-S,\\l, (28) 

quantifies the commute time in theory of random walks being equal to the 
expected number of steps required for a random walker starting ai i V 
to visit j € V and then to return back to i, 16J. 

It is important to mention that the cosine of an angle calculated in accor- 
dance to (|26p has the structure of Pearson's coefficient of linear correlations that 
reveals it's natural statistical interpretation. Correlation properties of flows of 
random walkers passing by different paths have been remained beyond the scope 
of previous studies devoted to complex networks and random walks on graphs. 
The notion of angle between any two nodes of the graph arises naturally as soon 
as we become interested in the strength and direction of a linear relationship 
between two random variables, the flows of random walks moving through them. 
If the cosine of an angle (|26p is 1 (zero angles), there is an increasing linear rela- 
tionship between the flows of random walks through both nodes. Otherwise, if 
it is close to -1 (tt angle), there is a decreasing linear relationship. The correla- 
tion is {tt/2 angle) if the variables are linearly independent. It is important to 
mention that as usual the correlation between nodes does not necessary imply 
a direct causal relationship (an immediate connection) between them. 
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3 Discussion and Conclusion 



In his speech to the meeting in the House of Lords, October 28, 1943, Sir 
Winston ChurchiU had requested that the House of Commons bombed-out on 
the night of May 10, 1941, be rebuilt exactly as before. "We shape our buildings, 
and afterwards our buildings shape us, " he said. 

A belief in the influence of the built environment on humans was common 
in architectural and urban thinking for centuries. There is a tied connection 
between physical activity of humans, their mobility and the layout of buildings, 
roads, and other structures that physically define a community [50^. Spatial 
organization of a place has an extremely important effect on the way people 
move through spaces and meet other people by chance [H]. It has been also 
shown in '24' that space structure and its impact on movement are critical to 
the link between the built environment and its social functioning. In particular, 
the reduce of movement in a spatial pattern is crucial for the decline of new 
housing areas. It is well known that the urban layout effects on the spatial 
distribution of crime [22] and poverty [53]. 

Equilibrium flows in transport networks can be investigated by means of 
Markov's processes defined on the dual graph representations of the network. 
Sustained movement patterns are generated by a subset of automorphisms of 
the graph spanning the spatial network of a city. This process is naturally 
interpreted as random walks. Random walks assign absolute scores to all open 
spaces in the city accordingly to the quality of paths they provide to random 
walkers and embed city space syntax into Euclidean space. 
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